On 

o 



^ 



X 
£ 



Universal amplitude ratios 

in the two-dimensional g-state Potts model and percolation 

from quantum field theoryQ 



b,c 



G. Delfino a and J.L. Cardy 

a Laboratoire de Physique Theorique, Universite de Montpellier II 

PL E. Bataillon, 34095 Montpellier, France 

b Theoretical Physics, University of Oxford 

1 Keble Road, Oxford 0X1 3NP, United Kingdom 

c All Souls College, Oxford 

<N 

We consider the scaling limit of the two-dimensional q-state Potts model for q < 
4. We use the exact scattering theory proposed by Chim and Zamolodchikov to 
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determine the one and two-kink form factors of the energy, order and disorder 
operators in the model. Correlation functions and universal combinations of critical 
amplitudes are then computed within the two-kink approximation in the form factor 
approach. Very good agreement is found whenever comparison with exact results 
is possible. We finally consider the limit q — > 1 which is related to the isotropic 
percolation problem. Although this case presents a serious technical difficulty, we 
predict a value close to 74 for the ratio of the mean cluster size amplitudes above 
and below the percolation threshold. Previous estimates for this quantity range 
from 14 to 220. 
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1 Introduction 

Two-dimensional integrable models enjoy a very particular role in the framework of quan- 
tum field theory. As a matter of fact, they provide the only known examples of non-trivial, 
interacting relativistic quantum theories to be exactly solved. As is well known, this cir- 
cumstance has to be traced back to the existence in these models of non-trivial integrals 
of motions resulting in a strong simplification of the scattering theory and in the possi- 
bility of determining the exact S'-matrix through bootstrap techniques. A large number 
of integrable models has been discovered and solved in this way over the last two decades 

If the theoretical relevance of integrable field theories is obvious, surely their effective- 
ness for the accurate quantitative description of interesting physical systems is not less 
important. In this respect, two-dimensional statistical mechanics offers a natural testing 
ground. Many statistical mechanical models are known to be exactly solvable in two 
dimensions || and are then natural candidates for a description in terms of integrable 
field theory in the scaling limit in which the correlation length becomes much larger than 
the lattice spacing. Moreover, one of the remarkable results of the recent investigations 
of two-dimensional quantum field theory is that integrability can be a property of the 
scaling limit of statistical models even in cases where it has never been found on the 
lattice. The critical Ising model in a magnetic field provides the most known example of 
such a situational^]. 

It is clear that, for the purpose of the quantitative study of statistical systems through 
the methods of quantum field theory, only universal quantities should be considered, i.e. 
those quantities which do not depend on the particular microscopic realisation of the 
system and are determined instead by its global features (essentially, internal symmetries 
and dimensionality). Universal properties emerge when the system approaches a second 
order phase transition point. For a magnetic system in zero external field, the usual 
characterisation of criticality is in terms of the classical thermodynamic quantities: in the 
vicinity of the critical point the specific heat, spontaneous magnetisation, susceptibility 
and correlation length behave as 

C~(A ± /a)r- a , 
M ~B (-rf , 

e=^T-", (i.i) 



x It is interesting to notice that for this specific case a lattice model has been identified which is in 
the same universality class than the Ising model in a magnetic field and is integrable pdj|. 



where we denoted by r the reduced temperature, r = a(T — T c ), a > 0, and the labels 
± on the critical amplitudes refer to the critical point being approached from above or 
from below. While the critical exponents are characteristic of the critical point and in 
d = 2 are determined by conformal field theory (CFT) || ||, the critical amplitudes 
carry information about the renormalisation group trajectory along which the critical 
point is approached and their determination requires a study of the system away from 
criticality. Critical amplitudes depend on metric factors and are not themselves universal, 
but universal combinations of them can by constructed which characterise the scaling 
region around the critical point. The following universal amplitude ratios are usually 
considered in the literature R 



a + /a. , r + /r_ , £ +/£ - , 



R c = A + T + /B 2 , R+ = A l l d Zt ■ (1.2) 



The scale factor independence of Rc and Rf is a direct consequence of the scaling and 
hyperscaling relations 

a = 2- 2(3 -7, 2-a = du. (1.3) 

The critical amplitudes can be expressed as moments of correlation functions of the 
spin and energy operators in the off-critical theory. As a consequence, the computation of 
the amplitude ratios in the framework of integrable field theory requires bridging the gap 
between the scattering theory, in terms of which the solution of the model is given, and 
the off-shell dynamics. With respect to this problem, the so called form factor bootstrap 
is the method which proved so far as the most effective. In this approach, correlation 
functions are expressed as spectral sums over n-particle intermediate states. The operator 
matrix elements entering the decomposition (known as form factors) are subject to a set 
of monodromy and residue equations ||, ||] which have been solved exactly for many 
integrable models. These equations, however, are fixed by the ^-matrix alone and do 
not distinguish between different operators. Although the selection rules coming from 
internal symmetries together with some minimality assumptions are sufficient in some 
cases to identify the form factor solutions corresponding to specific operators, it turns 
out that in general more information about the operator space has to be injected in order 
to handle this problem [10, 11 1. 



The operators which appear in physical applications are the scaling operators. Since 
they are naturally defined in the scaling limit towards the ultraviolet fixed point, it is 
not surprising that the constraints for their identification in the form factor approach 
come from the high energy asymptotics of the matrix elements. Two such constraints 
have been identified which are crucial for the results obtained in this paper. The first one 
relates the asymptotic behaviour of form factors of the scaling operators to their scaling 



dimension [ ID| ; the second constraint takes the form of a cluster decomposition of the 



matrix elements when the momenta of a subset of particles become much larger than the 



others. It was argued in Ref. Jl2| that the latter property is related to the decoupling of 
the theory into holomorphic and anti-holomorphic sectors in the ultraviolet limit. 

This level of understanding proved so far sufficient for the determination of the two- 
particle form factors of physically relevant operators in integrable theories. From the 
general point of view, the computation of the two-particle matrix elements amounts to 
fixing the initial conditions of the bootstrap procedure. The determination of the form 
factors with more than two particles through the residue equations is a mathematical 
problem which can be straightforward or extremely difficult depending on the degree of 
complexity of the underlying scattering theory. Although the problem has been solved in 



several specific cases (see e.g. [28, |9|, |T3||), no general method of solution is known for the 
case when the scattering allows for the exchange of quantum numbers among particles 
with the same mass. 

If in principle this circumstance severely restricts the range of applicability of the 
form factor approach for the exact computation of correlation functions, in practice very 
accurate results can be obtained anyway in a large number of cases. In fact, it has 
been by now verified for several integrable models that the spectral series over form 
factors is characterised by a remarkably fast rate of convergence (see, among others, 
Refs. |L4], [L3|, [15], [10]]). A theoretical justification of this property was proposed in ||T5| 



relying on phase space considerations and some peculiarities of integrable dynamics. 
As a matter of fact, it turns out that zeroth moments of two-point correlators can be 
computed within a typical accuracy of order 1%, including in the spectral series the one 
and two-particle contributions only. At fixed level of approximation, the accuracy rapidly 
increases if higher moments are considered in which the contribution coming from the 
short distances (more sensitive to the exclusion from the sum of many-particle states) is 
suppressed. 

It is the purpose of this paper to use the form factor approach to compute the uni- 
versal amplitude ratios for the g-state Potts model (q < 4) and the isotropic percolation 
problem. Although the latter does not fit immediately in the standard terminology of 
thermal phase transitions, it can be formally related to the q — > 1 limit of the Potts 
model. Through this mapping, the thermodynamic quantities listed in ( |1 . 1| ) become re- 
lated to the mean number of clusters, the percolation probability, the mean cluster size 
and the pair connectivity, respectively. 

Our starting point will be the exact scattering theory for the g-state Potts model 



proposed by Chim and Zamolodchikov in Ref. [16|. Due to the fact that the fundamental 



constraint of permutation symmetry which characterises the model is very naturally 



imposed on the interaction of the kinks interpolating between the q degenerate vacua 
of the spontaneously broken phase, the scattering theory is formulated at T < T c . Form 
factors over kinks have not been previously considered in the literature. Although most of 
the formalism used to deal with ordinary particles goes through with minor adaptations, 
some interesting new features appear. This happens, in particular, when operators which 
are non-local with respect to the kinks (e.g. the magnetisation operator) are considered. 
It turns out that the way non-locality is implemented in the low-temperature formalism 
is markedly different from that characteristic of the unbroken phase (see []nj for a review 
of the latter). 



The amplitude ratios ( |1.2|) involve both high and low-temperature amplitudes. We 
will work all the time in the low-temperature phase and rely on duality to get the infor- 
mation about the unbroken phase. In particular, the correlators of the magnetisation at 
T > T c will be obtained computing those of the disorder operators at T < T c . An inter- 
esting problem, however, arises at this point. As the form factors of the magnetisation 
and disorder operators are determined quite independently from each other, the relative 
normalisation between the two operators remains unfixed. This needs to be determined if 
the two operators are to describe the same physical quantity in the two different phases. 
The obvious solution to this problem would be to compute the two-point correlators of 
the two operators and to match the coefficients of their short distance asymptotics. In 
practice this cannot be done even in the cases in which all the multi-particle form factors 
are known, simply because nobody knows how to resum the spectral series. Once again 
the solution is provided by the asymptotic properties of form factors and in particular by 
the cluster factorisation mentioned above: a two-kink magnetisation matrix element has 
to factorise into the product of two one-kink disorder form factors, and this requirement 
fixes the relative normalisation. 

We will compute correlation functions and their moments in the g-state Potts model 
within the two-kink approximation in the form factor approach. When comparison with 
exact results is possible, the fast convergence of the spectral series is confirmed and results 
with the typical accuracy mentioned above are obtained. A serious technical difficulty 
however arises in the computation of the two-kink form factor of the magnetisation op- 
erator. Here, the conspiracy of non-locality and non-diagonal scattering results in a new 
form of monodromy problem expressed by a functional equation that we are not able to 
solve for generic values of q. The problem can be overcome for q = 2,3,4 and the complete 
list of amplitude ratios ( |1.2|) will be presented for this cases. For percolation, however, 
we can give accurate results only for some universal ratios. For the others we propose 
an extrapolation in q which, although naive, leads to results in good agreement with the 
available lattice estimates (series enumerations, Monte Carlo). Quite particular is the 



case of the ratio of the cluster size amplitudes r + /r_. Here the existing estimates range 
from 14 to 220, and our extrapolated value ~ 74 represents a considerable improvement 
in accuracy. 

The layout of the paper is the following. In the next section we briefly review the 
description of the scaling limit of the g-state Potts model as a perturbed CFT and the 
associated scattering theory. The one and two-particle form factors of the energy, order 
and disorder operators are determined in section 3. The results given by the form factor 
approach in the two-kink approximation for the central charge and scaling dimension sum 
rules, and for the universal amplitude ratios in the Potts model, are presented in section 
4. We discuss percolation in section 5 before summarising our results and making some 
final remarks in section 6. 

2 Scattering theory of the scaling Potts model 



The g-state Potts model |L7|, |lg] is the generalisation of the Ising model defined by the 
lattice Hamiltonian 

H = - J ^2 s s(x), s (y) > (2- 1 ) 

0:2/) 

where the sum is over nearest neighbours and the site variable s (x) can assume q possible 
values (colours). Clearly, the model is invariant under the group S g of permutations of 
the colours. In the ferromagnetic case J > we are interested in, the states in which all 
the sites have the same colour minimise the energy and the system exhibits spontaneous 
magnetisation at sufficiently low temperatures. There exists a critical temperature T c 
above which the thermal fluctuations become dominant and the system is in a disordered 
phase. If we introduce the variables 

a a (x) = 5 s{x):a , a = 1,2, ...,q (2.2) 

constrained by the condition 

i 

£>„(*) = 0, (2.3) 

the expectation values (a a ) differ from zero only in the low-temperature phase and can 
be used as order parameters. 

After defining x = e J l T — 1, the partition function of the model can be written in the 
form 

Z = Tr s Y[(l + xS s(x)tS{y) ) . (2.4) 

0,s/) 



A graph Q on the lattice can be associated to each Potts configuration by drawing a 
bond between two sites with the same colour. In the above expression, a power of x is 
associated to each bond in the graph. Taking into account the summation over colours 
one arrives to the expansion [III 



Z = J>*V\ (2.5) 

Q 

where iV^ is the total number of bonds in the graph Q and N c is the number of connected 
components (clusters) in Q (each isolated site is also counted as a cluster). In terms of 
the partition function ( |2.5| ) the g-state Potts model is well defined even for noninteger 
values of q. 

In two dimensions, the Potts model is known to undergo a first order phase transition 
at T = T c for q > 4 M. For q < 4, however, the transition is continuous and the critical 
point can be described by a CFT. Relying on the knowledge of some critical exponents 
20H , Dotsenko and Fateev identified the central charge of this CFT to be [^1 



(2.6) 



t(i + l)' 

where the parameter t is related to q by the formula 

^2 8ill ^i>. (2.7) 

The scaling dimension x a of the magnetisation operators cr a {x) (the continuous version 

of the quantities ( |2.2|) ) is identified with that of the primary operator (fi(t-i)/2,(t+i)/2 in 

the CFT 

= («-!)(« + 3 ) 

8t(t + l) V ; 

The energy density operator e(x) (~ ^2 y $s(x),s(y) on ^ ne lattice) coincides with the pri- 
mary operator 02, i with scaling dimension 

*. = \ (l + |) . (2.9) 

In view of these identifications, the continuum limit of the g-state Potts model (q < 4) 
is described by the perturbed CFT 

A = A C ft + t d 2 xe(x), (2.10) 

where Aqft denotes the critical point action. Since a series of non-trivial integrals of 
motions is known to survive the deformation of a CFT by the operator 02,i [@], the 



off-critical theory ( p.lOj ) is integrable. This circumstance was exploited in Ref. [|TI| to 
propose an exact scattering theory which we now briefly review 2 . 

The low-temperature phase of the model is characterised by the presence of q degen- 
erate vacua that we label by the index a = 1,2, ... ,q. The elementary excitations are 
then provided by kinksQ K a p{9) interpolating between the two vacua a and (3 (a ^ (3). 
The space of physical asymptotic states consists of multi-kink configurations of the type 
K aoai {9i)K aia2 {9 2 ) ■ ■ ■ K an _ lCtn (Q n ) («j ^ ctj+i) interpolating between the vacua «o and 
a n . As a consequence of the invariance under permutations, all the n-kink states fall 
into two topological sectors: the "neutral" sector, corresponding to a = a n , and the 
"charged" sector, corresponding to ao ^ a„. 

Integrability implies that the scattering processes are completely elastic and factorised 
into the product of two-kink interactions. Since topological charge is conserved, an 
outgoing two-kink state can only differ from the ingoing one by the vacuum state between 
the kinks. Hence, the two-kink scattering can formally be described through the Faddeev- 
Zamolodchikov commutation relation 

K^e^K^i) = J2 si s f3 (e 12 )K a s(e 2 )Ks P (e 1 ) , (2.11) 

S^a,f3 

where 6*12 = Q\ — 9 2 and 5^ (#12) denotes the two-body scattering amplitude (Fig. la). 
.Sg-invariance reduce to four the number of independent amplitudes, two for the charged 
and two for the neutral topological sector 

K a7 (9 l )K 70 (9 2 ) = S Q {9 l2 )Y, K as{02)K s ^9 x ) + S l (9 12 )K a7 (9 2 )K 7/3 (9 1 ) a + (3 
K ai (9 x )K ia (9 2 ) = 32(9 13 )^2K aS (e 2 )Ks a (e 1 ) + S-^K^^K^) , (2.12) 

<¥7 



Using the commutation relation Q2.11|) twice one obtains the unitarity constraint 



e 

which amounts to the set of equations 

(q - 3)S (9)S (-9) + S 1 (9)S 1 (-9) = 1 , (2.14) 

2 Notice that a different 5-matrix for the 02, 1 perturbation of minimal models was determined by 
Smirnov |22] . We use here the scattering description of Ref. ]16| because it is more suitable for analytic 
continuation in q. Both descriptions must lead to the same results for the correlation functions. See 123] 
for a detailed discussion of this point in the case of the 4>i,3 perturbation. 

3 We parameterise on-shell momenta as p M = (mcosh^msinhfl), where m denotes the mass of the 
kinks. 



(q - 4)S (6)S (-6) + So(6)S 1 (-9) + S 1 (9)S (-9) = , (2.15) 

(q - 2)s 2 {e)s 2 (-e) + s 3 (9)s 3 (-9) = 1 , (2.16) 

(q - 3)S 2 {9)S 2 {-9) + S 3 (9)S 2 (-9) + S 2 (9)S 3 (-9) = . (2.17) 

Crossing symmetry provides the relations 

S (9) = S (in - 9) (2.18) 

S t (9) = S 2 {m - 9) (2.19) 

S 3 {9) = S 3 (m - 9) . (2.20) 

Using these constraints together with the Yang-Baxter and bootstrap equations (that 
we do not need to reproduce here) the following expressions for the four elementary 
amplitudes were determined in Ref. [M 



S (9) = sinh Afl sinh A(fl^) / a*n 

V; sinh A (9 - 2f) sinh A (0 - f ) V™7 



sin 2|* sinh A(0 - ivr) /A0\ 
sin ^ sinh A (0 - ^f ) \Jk ) 



sm = :" j : ?; S n; (2.22) 



sin^ sinhA# /A0\ 

Safe) = j— 2 7 -r n — 2.23 

V ; sin^ sinhA(0-f) \iir J V ! 



sin A7r ^ ( \9\ 

3 

where A is related to g as 



s 3 (0) = -^rn — , (2.24) 

sm ^ V m J 



Vg = 2 sin ^- , (2.25) 



and 



r(l-x)r(l-A + x)r(|A-x)r(|A + x) -Fr , , , 

n(x) = r4 Mr r— tI c TTn*(x)n fc (A 



x 



n fc (x) 



T(l + x)T(l + X-x)T (|A + x) r (|A - x) A ., 

T(l + 2A;A - x)r(2A;A - x)T [l + (2k - ±) A - x] V [(2k + f ) A - x] 
T[l + (2fc + 1)A - x]r[(2fc + 1)A - x]T [1 + (2k - §) A - x] f [(2fc + f) A - x] 

We quote here the following integral representation of the function Il(x) which will be 
useful in the following 

nffl^"M^", (2.26) 

\iit J sinh X(9 — in) 

r 00 dx sinh f (1 - i) - sinh f (i - |) ,xfl / N 

4(0) = / 2 V . ¥n rf^ — sinh — . 2.27 

Jo x smh ^- cosh | 27r 

It is easily seen that the function Tl(\9/m) is free of poles in the physical strip Im# e (0, n) 
for q < 3 (i.e. A < 1). Hence, in this range of q the only poles of the scattering amplitudes 



in the physical strip are those located at 9 — 1m j3 and 9 = in/3 and correspond to the 
appearance of the elementary kink itself as a bound state in the direct and crossed 
channel, respectively. The residues 

Res0=2i7r/3»S'o(#) = —Rese =in / 3 S (9) = 



Res, 



Uw/3 
■e=2iw/3 



Sl (9) = -Res $=i7T/3 S 2 (9) = i(r« K y, 



determine the coupling at the three-kink vertex (Fig. 2a) 



KK 



1 2ttA 
A Sin ^ 



e ^f) 



1/2 



(2.28) 



For q > 3 (A > 1) a direct channel (positive residue) pole located at 9 — 2in, 



2 V A 



(2.29) 



enters the physical strip in the amplitudes S2(9) and S 3 (9). Such a pole must be accord- 
ingly associated to a (topologically neutral) kink-antikink bound state B with mass 



rriB = 2m cos k . 



(2.30) 



Of course, the amplitudes S\(9) and Ss(9) exhibit the corresponding crossed channel 
(negative residue) pole at 9 — m — Iik. The coupling at the kink-kink-bound state 
vertex (Fig. 2b) is given by 



Fkk = [-i^s e= 2 iK S 3 (9)} 



1/2 



1 47rA sin Xn 



— sin r— 

A 3 sin 4? 



,A(2in) 



1/2 



(2.31) 



The amplitudes Skb{9) and Sbb(9) (Fig. lb,c) describing the kink-bound state scattering 
and the bound state self-interaction are determined by the bootstrap equations 

Skb{9) = {q- 2)S 2 {9 - ik)S x {9 + in) + S 3 (9 - ik)S 3 {9 + in) , 

Sbb(9) = S BK (9 - ik)S bk {9 + in) , (2.32) 

and read 



Sbk{9) = ti- K / 7T (9)t2/3~K/Tr(9) , 
Sbb(9) = t2/ 3 (9)t 1 _ 2K / 7r (9)t 2 /3~2K/Tr(9) , 



in terms of the building blocks 



tanh h(9 + in a 

t a (9)- 



tanh g (9 — ma) 



(2.33) 



(2.34) 



The poles located at 6 = i(jt — k) in Sbk and at 6 = 2n/3 in Sbb are bound state poles 
corresponding to K and B, respectively. The coupling at the BBB vertex (Fig. 2c) is 



1/2 



^BB = \—i^Se=2i-K/2,SBB{0)\ 

1/2 

.(2.35) 



2^3 cot - ( 1 ) cot 7T ( — ) tan - ( 4 - - ) tan - ( 5 - - 

2 V XJ V2A3/ 6 V A 7 6\ A 



It can be shown [M that, at least in the region q < 4 we are interested in, the remaining 
poles in the amplitudes Skb and Sbb are associated to multi-scattering processes rather 
than to new particles^. Hence, the elementary kinks and their neutral bound state B are 
the only particles entering the spectrum of the theory in this range of the parameter q. 

3 Form factors 

The two-kink form factor of an operator $(x) (Fig. 3) is defined as the matrix element^ 

0*12) = K\®mK ai {6x)K lP {6 2 )) , (3.1) 

where |0 a ) denotes the vacuum state in which all the sites have colour a. The fundamental 
equations constraining this matrix element come from the requirements of unitarity and 
crossing symmetry. The unitarity equation for form factors follows immediately from 

Eg- (TO 

*£»(*) = E S l» F ^{-0) ■ (3-2) 

The crossing equations read (see the Appendix) 

F* ia (6 + 2tn) = K% 7 (-6) , (3.3) 

for the neutral sector, and 

F% p {B + 2m) = F% p {-B) , a ? (3 (3.4) 

for the charged sector. Let us consider how Eqs. ( |3.2|) -( J3T4"1) specialise for the different 
operators we are interested in. 

4 Multi-scattering singularities usually show up as higher order poles in two dimensions (they lead 
to anomalous thresholds in four dimensions) p5[ . In the present case a simultaneous vanishing of the 
residues reduces the singularities to simple poles |24| . 

5 We will consider only operators which are scalar under Lorentz transformations. Their matrix 
elements depend on rapidity differences only. 
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3.1 Energy operator 

The energy operator e(x) is the operator which perturbs conformal invariance in the 
action ( |2.10| ). It is then proportional to the trace of the stress-energy tensor Q(x) 



G{x) = 2tt{2-x £ )ts(x) 



(3.5) 



In view of this relation, we will mainly refer to Q(x) in the following. 

Being invariant under the S q symmetry, 0(x) couples to the neutral sector of the 
space of states. Moreover, its two-kink form factors 



FV(6 12 ) = (0 a \Q(0)\Kcrr(0i)K ia (9 2 )) 



(3.6) 



do not depend on the choice of a and 7. The unitarity and crossing equations can then 
be written in the form 



where 



Since 



F B (6) = A(6)F°(-6) 
F @ (6 + 2nt) = F @ (-t 



A(e) = {q -2)s 2{ e ) + sm = si ^ + m l m 



sinh X(9 — in) 



£{6) = exp 



9e{x) 



00 dx . . , x6 
g £ {x) sinh 



uo 



x 



171 



sinhfcosh(|-l)f 
sinh ^- cosh | 



sinh A (8 + in a) 
sinh A(^ — ina) 



exp 



■g a (x)smh 



lyo 



x 



ITT 



9a(x) = 2 
it is easily checked that the function 



sinh (A ~ a ) 
sinh^ 



x 



. . ,9 ( f°° dx g x (x) + g £ {x) . 2 {m - 9)x 
^a{9) = —1 smh - exp < / — sin 



x sinh x 



2tx 



(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 
(3.12) 



, A < 1 (3.13) 



solves the system of functional equations (|3.7| ). We note here for later convenience the 
asymptotic behaviour 



F A (8)~ exp (1-^10. 



-00 . 



(3.14) 
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The integral in (|3.13| ) is convergent on the real #-axis for A < 1. For A > 1 one needs to 
use the analytically continued expression 

^f^^F^W^' A>1 (3 - 15) 

t/ m . . , 9 J f°° dx gi-i/x(x) + g E (x) . 2 (in-e)x\ 

T(0) = -*sinh-expjy o — sm —^—j , (3-16) 

which explicitly exhibits the pole at 9 = lin in the physical strip corresponding to the 
kink-antikink bound state B (Fig. 4a). Taking into account the normalisation condition 

F e (nr) = 27rm 2 , (3.17) 

we can write 

F e {6) = 27im 2 F A {6) , (3.18) 



. T B O-AK V I s[n2K T B 



F% = (0 tt |6(0)|B> = -^^Res e=2 i^(e) = —f- -^ . (3.19) 



Let us finally determine the matrix element F BB {0i 2 ) = (0 a \Q(0)\B(9 1 )B(9 2 )) . It has to 
satisfy the monodromy equations 

F® B (9) = S BB {9)F® B {-9) , 

F% B (6 + 2m) = F§ B (-e) . (3.20) 



The solution to these equations with the expected pole structure is 


2 

where 



FbbW) = [acosh9 + b]smh-K 2/3 {9)K 1 _ 2K/lT {9)K 2/3 „ 2K/ M , (3.21) 



, f°° dx cosh (a- \)x 9 (w - 60x | 

^W = — r^ exp 2/ \ , y sin 2 l ; . 3.22 

cosh y — cos 7ra | J x cosh | smh x In 

The coefficients a and 6 in ( |3.21| ) are uniquely determined by the residue equation 

- i Res e=2W3 i 71 I B (^) = r| B F| , (3.23) 

and the normalisation condition 

F® s (ivr) = 2vrm| . (3.24) 

Having computed F B and F BB (9), we can determine the vacuum expectation value F® = 



(0 Q |G|0 Q ) through the asymptotic factorisation relation |T2 



lim F* B (9) = ^ • (3.25) 
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We find 

n tc sin 7JT n , 

in perfect agreement with the result obtained in J26] using the thermodynamic Bethe 
ansatz (TBA). 

3.2 Disorder operators 

The theory contains disorder operators which are dual to the q — 1 independent order 
operators cr a (x) and have the same scaling dimension ( |2.8| ). Following the usual inter- 
pretation of the disorder operator as the kink creation operator, we denote by /ig(x) the 
disorder operator which, acting on the vacuum state |0 Q ), creates states interpolating 
between the vacua a and (3 (/3 ^ a). The two- kink form factor 

(0 a \^(0)\K a ^9 1 )K lS (e 2 )) = 5p S F»(0 12 ) , (3 ^ a (3.27) 

does not depend on the intermediate index 7 and satisfies the unitarity and crossing 
equations 

F fi (9) = Z{9)F»(-9) , 

F"(0 + 2m) = F"(-0) , (3.28) 

with 

sinh A (9 + — ) 
E(0) = (g - 3)S O (0) + &(0) = sinhA L_gl *(*) • (3-29) 

Hence, it can be written as 

F " {6) = Z " sinh \ (9 + f^sinh § (0 - f) ' (3 - 30) 

where 

pm . . ,0 / f°° dx g 2/3 (x) + g £ (x) . 2 (in-6)x \ 

Fx{V) = -1 sinh.- exp< £ — sin } , (3.31) 

2 [Jq x smni Itx J 

is the solution of the system ( |3.28| ) without poles in the physical strip; the pole at 

9 = 2iir/3 corresponding to the </2 3 -property of the kinks has been inserted explicitly. 

The function F%(9) behaves asymptotically as 

F s (0) ~ exp M - ^ J 9 , — > +00 . (3.32) 

The one-kink form factor 

(0 a \^(0)\K ay ) = 8p 7 F£ , (3 ± a (3.33) 

is given by (Fig. 4b) 

K = ^Res e=2W3 i^W = - 4F fZ /3) Z » ■ ( 3 - 34 ) 

11 KK VdL KK 
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3.3 Magnetisation operators 

The magnetisation operators a a (x) couple to neutral states. Since a a carries an index, its 
matrix elements are not in general invariant under permutations. In the two-kink form 
factor we distinguish the three components 

^ a \a (Q)\K ai (9 x )K ia (9 2 )) = 8 0a F[(9 12 ) + 5 Pl F°(9 l2 ) + (1 - S Pa - 5^)F^(9 12 ) , (3-35) 

which eq. fl2.3|) relates as 

Fi(P) + F? (9) + (q- 2)F°(9) = . (3.36) 

Eq. (|3.2|) provides the relations 



F?(9) = A(9)F?(-9), 

F°(9) = S 3 (9)F 2 (-9) + (q- 2)S 2 (9)F°(-i 

F°(9) = S 2 (9)F°(-9) + [(q - 3)S 2 (9) + S 3 



Fs(-t 



(3.37) 



The crossing equations 



F?(9 + 2tn) = F°(-9), 
F°(9 + 2tn) = F°(-9), 



(3.38) 



are an immediate consequence of (|3.3|). Recalling that the function F A (9) defined in 
( |3.13|) satisfies Eqs. (|3.7|), we parameterise F°(9) as 



F[{9) = £l(9)F A (9) ; (3.39) 

then, the first of (|3.37| ) is automatically fulfilled provided the function £1(9) satisfies 

£1(9) = £l(-9) . (3.40) 



Simple manipulations involving Eqs. ( |3.36| )-( j338| ) provide the basic equation that £1(9) 
should satisfy together with the previous one: 



£1(9) 



S 3 (6) 
S 2 (9) 



- 1 



£l(2m + 



A(9) 



£l(2m - 



or, more explicitly, 

£1(9) -- 



sinh X(9 — in) £l(2in 



S 2 (9) 
+ sinh X(9 + in) £l(2in 



2 cos ^ sinh X9 



(3.41) 



(3.42) 



Unfortunately, we do not know how to solve this equation for generic values of A. Before 
turning to the solution for specific values of this parameter, we collect some additional 
physical information about the magnetisation matrix elements. 
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The order operators cr a (x) are non-local with respect to the disorder operators t/,p(x) 
which interpolate the kinks. This kind of mutual non-locality is known to lead to the 



presence of an "annihilation pole" at 9 = m in the two-particle form factor || [14], ^7 
We show in the Appendix that in our case the residue on such a pole takes the form 

-iRes 01 _ e2=i7r (O a |a 7 (O)|^ a/3 (^)^ a (0 2 )) = (0> 7 |0 Q ) - (O^a^) . (3.43) 

Denoting 

F^ = (0 a \a a \0 a ) , (3.44) 

and taking into account the constraint fl2.3|), the order parameter can be written as 



(0aN0«) = -^ r (^ 7 a-l) • (3.45) 

q — 1 
Hence, for the different components of the two-kink form factor we conclude 

-iRes e=i „Ff{9) = -iRes 9=iw FZ(9) = -37^ 

-iRes e=in F£(9) = (3.46) 

For A > 1 the magnetisation has a one-particle form factor on the bound state B. 
Denoting 

F% = (O a \(r a (p)\B), (3.47) 

we can write 

(0> 7 (0)|5> = -^- (q5 ia - 1) . (3.48) 

q — 1 

These one-particle form factors can be obtained from the two-kink matrix elements 

through the residue equation 

-iRes ei -e^2i K ^a\^)\Kap{ex)Kp a {92)) = I** (O a |<r 7 (0)|B) . (3.49) 

The last two equations, together with the crossing relation ( p.38 ), immediately lead to 
the following identity for the function Q(9) introduced in (|3.39|) 

VL(2m - 2m) = ft(2m) . (3.50) 

q-1 

Although obtained for A > 1, this relation is expected to hold for generic values of the 
parameter. 

The matrix elements of scaling operators in unitary theories must satisfy the asymp- 



totic bound obtained in Ref. p0| . For the kink- kink form factor of the magnetisation 
operators this reads 

lim F 1 CT (9) < constant e x ° e/2 . (3.51) 

9 — >+oo 
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The functional equation (|3.42|) becomes trivial at the points q = 2,3,4. Let us 
investigate these cases in more detail. 
q=2 

For this value of q (corresponding to A = 3/4) the solution to Eqs. (|3.40p and ( p.42 ) 
satisfying the asymptotic bound ( |3.51| ) is simply 



Q(9) 



cosh 7 



(3.52) 



where the normalisation has been fixed using Eq. (|3.46|) . The identity (|3.50| ) is automat- 
ically satisfied. Since F\(6) = — 2sinh#/2, one obtains 



F?(8) = -F°{0) = iF ff tanh 



6 



(3.53) 



which is the well known result usually obtained in the high temperature formalism [28, 14 



q=3 

For A = 1 Eq. ( |3.42j ) reduces to 



ft(0) = _[fi(0 + 2m) + Q(9 - 2in)} 
Together with Eqs. (|3Tl0| ), ( p3I]) and (|3T46| ), it fixesQ 

n(e)-- 



y/3 cosh ■ 

•— — Fn 



cosh ■ 



Again Eq. (|3.50|) is automatically satisfied. 



(3.54) 



(3.55) 



q=4 

For A = 3/2 Eq. (ggg) becomes 



£1(6 + 2iit) = £1(6 - 2m) . 



(3.56) 



This time Eq. ( |3.50|) has to be enforced together with the other constraints in order to 
fix the solution 



0.(6) 



2F CT cosh | + y/3 



(3.57) 



3^3 cosh f 

We conclude this section considering the problem of the relative normalisation of the 

order and disorder operators we discussed in the introduction. Solving Eq. Q3.42 ) for 

6 — > oo one easily finds that fl(6) has to behave in this asymptotic limit as e ve , with 

6 An equivalent result can be obtained working with the high temperature scattering theory which is 



defined in terms of particles rather than kinks |29, 30 
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r] = A/3 + k or r] = 2A/3 + k {k integer). The solutions obtained at q = 2, 3, 4 indicate 
that the correct choice is 77 = 2 A/3 — 1. Combining this result with the asymptotic 
behaviour ( |3.14|) for F\(8), we conclude that F°(8) goes to a constant in the limit of 



large rapidity difference. Hence, according to the discussion of Ref. |12]1 , we expect the 
following asymptotic factorisation equation to hold 



lim |Ff 






(3.58) 



In the next section, this relation will enable us to express the critical amplitudes in terms 
of a single arbitrary normalisation constant (say Fq) which in turn cancels out when the 
universal amplitude ratios are considered. 



4 Correlation functions and amplitude ratios 

The knowledge of the form factors allows to express correlation functions as spectral sums 
over intermediate asymptotic states. For the two-point euclidean correlator of two scalar 
operators $i(x) and $2(2;) we have 



^ Je 1 >...>e n 2tt 2tt 



(0)|n>(n|$ 2 (0)|0)e- 



-\x\En 



(4.1) 



where E n denotes the total energy of the n-particle state \n). 

Clearly, the above expression is a large distance expansion: while the intermediate 
states with the lowest total mass provide the dominant contribution for large separations, 
in principle the whole series should be resummed in order to reproduce the correct ul- 
traviolet behaviour. Relying on the properties of fast convergence of the spectral series 
we mentioned in the introduction, we consider partial sums of the series truncated at the 
level of the two-kink intermediate state (the relevant form factors have been computed 
in the previous section). For example, the expansions for the two-point correlators of the 
order and disorder parameters read 



(0 a Mx)a 7 (0)|0 a > 



(qSpa- l)(qS. 



7 tt 



1: 



I P°"l 2 



1 

VT ./ 



(q-iy 

dd f^(29) K (2m\x\ cosh 



H(q-3)^K (m B \x\) 



T\ 



+ . 



+ 

(4.2) 



(0 Q | / u /3 (x)/i 7 (0)|0 a ) 



+ 



8fa 
Q- 



l£_fc I 

71 

2 rco 



7i z 



K (m\x\) + 
d9\F> 1 (2e)\ 2 K (2m\x\coshl 







+ 



/3, 7 ^o(4.3) 
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where we introduced the function 

/#r(0) = [Mt«(? - !) - W 1 " M " (1 - <W)<Ml^iW + 

+ (1 - 6fr)(l - 6 0a )[F 2 °(9)F°(-0) + FZ(6)FZ(-6) + {q- 3)|F 3 *(0)| 2 ] ,(4.4) 

and the step function H(y) which equals 1 for y > and is otherwise. The terms we 
omitted in the r.h.s. are of order e -3m l x l for |x| large (at least for q < 3). We show 
immediately that this level of approximation is sufficient to obtain remarkably accurate 
numerical results for the quantities we need to evaluate in this paper, namely moments of 
two-point correlators of the type J d 2 x|x| p ($i(x)$ 2 (0)) c . An interesting check is provided 



by the following sum rules for the central charge of the ultraviolet CFT [31| and the scaling 
dimension of the magnetisation operator |T2] 



c = ^ J d 2 x\x\ 2 (o a \e(x)e(o)\o a )c 

x ° = ~ o /n i in \ [d 2 x{O a \Q{x)a 7 {0)\O a } c , (4.5) 

27r(U Q |cr 7 |U a ) J 

where (• • -) c denotes connected correlators. The result obtained for c as a function of q 
using the truncated spectral expansion of the trace-trace correlator is shown in Fig. 5 and 
compared with the exact value 02.6D I 7 !- The numerical results obtained for c and A CT for 
q = 2,3,4 are listed in Table 1. Since the energy operator in the two-dimensional Ising 
model couples to the two-kink state only, the results obtained for q = 2 are exact. In the 
other cases the observed deviation from the exact result is of few percent at most. 

The critical amplitudes ( |1 . 1| ) are linked to the off-critical correlators by the relations 



C= J d 2 x(O a \e(x)e(0)\O a ) c , (4.6) 

M = (0 Q |a Q |0 Q ), (4.7) 

X= I d 2 x(0 a \a a (x)a a (0)\0 a ) c , (4.8) 



1 J d 2 x \x\ 2 (0 a \a a (x)a a (0)\0 a ) c 



4 Jd^x(O a \a a (x)a a (0)\O a ) c ' V ' ; 



7 It can be appreciated from the figure that the values of c obtained from this first contribution to 
the spectral sum are slightly larger than the exact result in the range 1 < q < 2. At first sight this 
may seem strange since each intermediate state which remains to be included into the sum will give 
a positive contribution (an integral over the modulus square of the corresponding trace form factor). 
Notice however that the number of three-kink intermediate states in the neutral topological sector is 
(q — l)(q — 2) and becomes negative when f < q < 2! 



The dimensional parameter entering our S-matrix approach is the mass m of the kink. 
It is related to the reduced temperature r appearing in ( |1.1| ) as^ 

m = rriQT u . (4-10) 

Equation (|4.9| ) defines the "second moment" correlation length. In the literature the so- 
called "true" correlation length £ t is often considered which is defined through the large 
distance asymptotic decay of the spin-spin correlators 

(a(x)cr(O)) ~ exp(-\x\/£ t ) , M — > oo . (4.11) 



It follows from fl£|) together with duality that & = l/m at T > T c . At T < T c , Eq. (gj) 
implies instead £t = l/2ra for q < 3, and £ 4 = l/m^ for 3 < q < 4. 

Since the zeroth moment of the energy-energy correlator appears in the scaling di- 
mension sum rule 

%e = - * , n , [d 2 x(0 a \Q(x)e(0)\0 a )c, (4.12) 

27r(0 Q |£|0 Q ) y 

the specific heat amplitudes A± can be computed exactly as[] 

A± = -a(l - a)(2 - a) (— Y ^- , (4.13) 

\m J An 

where F® is given in ( p.26[ ). The equality of A + and A_ follows from the fact that the 
energy operator e(x) simply changes sign under the duality transformation exchanging 
the low and high temperature phases. 

Remembering the definition ([3.45 ), the magnetisation amplitude B can be written as 

B=H*>». (4.14) 

The susceptibility and correlation length amplitudes can be evaluated using the expan- 
sions ( |4.2| ) and ( |4.3| ) for the correlators and the corresponding form factors. By duality, 
the high temperature amplitudes are obtained substituting o a by /ip {(3 ^ a) in Eqs. 
(4.8) and (|4.9| ). Some of the values obtained in this way for integer q are listed in Table 2. 



In particular, the results for £q mo show that the "second moment" and "true" correlation 
lengths differ very slightly from each other at T > T c . 

Table 3 collects the results we find for the amplitude ratios at q — 2,3,4. Excepting 
A+/A-, we are not aware of previous reliable estimates of these ratios for the cases 

8 The dimensionless number too is exactly computable through the TBA [£6|; we will not need it here. 
9 The specific heat diverges logarithmically in the Ising model (a — 0) and the definition of the 



result A± — to,q/27t 



amplitudes is accordingly modified to C ~ — A± lnr. In the limit q — ► 2, Eq. ( 4.13 ) gives the correct 
ll 
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q = 3,4. At q = 2, however, comparison with the known exact results |32|, Q (see also 



T + /r_ = 37.69365.. , 
R c = 0.318569.. , 

£ + /^=3.16.., (4.15) 

further confirms the remarkable accuracy of the results yielded by our two-kink approxi- 
mation. 



5 Percolation 

Percolation is the purely geometrical problem (no temperature involved) in which bondsQ 
are randomly distributed on a lattice with occupation probability p |35| . A set of bonds 
forming a connected path on the lattice is called a cluster. There exist a critical value 
p c of the occupation probability above which an infinite cluster appears in the system; 
p c is called the percolation threshold. If N is the total number of bonds in the lattice, 
the probability of a configuration with iV& occupied bonds is p Nb {l — p) N ~ Nb . Hence, the 
average of a quantity X over all configurations Q is 

(X) = J2Xp Nb (l-p) N - Nb . (5.1) 

Q 

Let x and y denote the positions of two bonds on the lattice and consider the function 
C(x, y) which takes value 1 if x and y belong to the same cluster, and otherwise. Then, 
the function g(x, y) = (C(x, y)) is the probability that x and y belong to the same cluster 
and is called pair connectivity. If we denote by P the probability that a bond belongs to 
the infinite cluster (P = at p < p c ), clearly we have 

lim g(x,y) = P 2 . (5.2) 

a; — y| >oo 

Since J2 V C(x, y) counts the total number of bonds in the cluster x belongs to, the mean 
cluster size can be obtained as 

S = ^2g e (o,y), (5.3) 

y 
where the subscript c means that the connected part of the pair connectivity is taken in 
order to get rid of the contribution coming from the infinite cluster at p > p c . The second 



10 We refer here to bond percolation; in site percolation, sites rather than bonds are occupied with 
probability p. Of course, all universal results are independent of this distinction. 
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moment correlation length is also naturally denned in terms of the pair connectivity as 

? 2d J2 x 9c(x,0) ■ l • ' 

The following relation with the g-state Potts model is particularly important for the 
theoretical study of percolation processes [ S5| . Remembering ( |2.5|) , the average of a 



quantity X in the g-state Potts model can be written as 

{X) q = Z- x Y,X<F c x Nh \ (5-5) 



Q 

hence, it is sufficient to make the formal identification x = p/(l — p) to see that (X)i 
coincides with the percolation average ( |5.1| ). For example, the mean cluster number in 
percolation can be expressed as 

(N c ) = lim — — = lim . (5.6) 

9—1 dq 9—1 q - 1 

To proceed further with this mapping we need a representation of pair connectivity 
in the Potts model formalism. For this purpose observe that the insertion of a delta 
function S s ( x ) a in a Potts configuration fixes to a the colour of the cluster x belongs to. 
Hence, at T > T c , 

\ys(x)alq — ) 

q 

(6 a (x)a$s(y)a)q = ~ 9g(x, y) + — [1 - g q (x, y)\ , (5.7) 

where g q (x,y) is the probability that x and y belong to the same cluster. Using the 
definition (|2.2|) we immediately find 

q-l 

{a a {x)a a {y)) q = — — g q (x, y) , (5.8) 

from which we see that the pair connectivity in percolation can be obtained as 

g(x,y) = lim (a a (x)a a (y)) q . (5.9) 

q ►! q — 1 



Comparison with ( |5.2|) , ( |5.3| ) and ( |5.4j) provides the following relations with the magneti- 
sation, susceptibility and correlation length in the Potts model 

M 



P = lim 



S = lim , 

9—i q-l 

£=£| ? =i- (5-10) 
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Together with (|5.6|) , they imply the following critical behaviour near the percolation 
threshold 



V*c, 


- ~ A± 


:\Pe- 


n\ 2 ~ a 


N 


P\ j 


p ~ 


B(p- 


■Pcf 


i 


S ~ 


f±|Pc 


-p\~ 


7 

5 



i^&bc-pr*, (5.11) 

where the critical exponents are those of the Potts model evaluated at q = 1, and the 
amplitudes are related to the Potts amplitudes asQ. 

7 i- -< 1 ±'0 
A+ = lim r ; — , 

9— »i(g-l)a(l-a)(2-a) 



.B = lim 



fir, 



./3 







9— »i W 
lim 



iW 



) 



9— >i g - 1 

&=^o _ 1 g=1 ; (5-12) 

here we introduced the non- universal positive constant r entering the relation r ~ T (p c — 
p) . It follows from these relations that the following combinations of critical amplitudes 
in percolation 

A + /A_ = lim A + /A_ , 

q >1 

f + /f_= iimr + /r_, 

q >1 

SO /SO = SO /SO 1 9=1 J 

2 _ l;„ ^ 



Rc = a(l - a) (2 - a)A+r + / J B 2 = lim 



4 = [a(l - a) (2 - a)A,.] 1/d & = Hm _ * , (5.13) 

are universal and can be computed from the q — > 1 limit of the Potts amplitude ratios. 
Let us see which results we can obtain for percolation in d = 2 from our study of the 
g-state Potts model of the previous sections. From (|4.13|) we find 



2^-2-a 



.4 



± 



m 5 r o 



2v/3 



(5.14) 

9=1 



The negative sign of this amplitude agrees with the series result of Domb and Pearce 
[ J3?H which is listed in Table 4 together with other series and Monte Carlo estimates of 



11 Notice that the labels + and — refer to p < p c and p > p c , respectively. 
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percolation amplitudes^. The equality of A + and A_ follows from duality which is also 
a crucial ingredient of the Domb and Pearce lattice calculation^. When combined with 
the value of £q~ quoted in Table 2, (|5.14|) gives the result (a = —2/3 at q — 1) 



RJ - 0.926, (5.15) 

which we think substantially improves the value around 1.1 one can extract from the 
lattice amplitudes in Table 4. 

We cannot compute directly the ratios involving the amplitudes B, T_ and £q simply 
because we are not able to solve the functional equation ( |3.42| ) for q = 1. What we can 
do is to attempt a naive quadratic extrapolation at q = 1 of the results obtained for 
q = 2,3,4. Since all the results following from the scattering theory are analytic in A, 
we perform a quadratic extrapolation in this variable^. From the values of Table 2 one 
extrapolates T + ml/B 2 ps 3.5 at q = 1, which in turn leads to Rc = 4:0T + ml/27\^3B 2 ps 
3.0. Similarly, the results of Table 3 lead to &/& « 3.76 and f + /f_ ps 74.2. The 
extrapolated values for £q/£q and Re compare quite well with the Monte Carlo result 
£+/£- = 4.0 ± 0.5 of Ref. ||] and with the estimate R c « 2.7 - 2.8 we deduce from 
Table 40. 

The status of the lattice estimates of the ratio r + /T_ is extremely controversial. While 
there exists a substantial agreement (within 20-30%) on the value of T + , the amplitude 
T_ turns out to be very small and difficult to determine. This resulted in a series of 
estimates for the ratio ranging from 14 to 220 (see 0)! When a value around 200 seemed 
to be accepted, the authors of Ref. |39) obtained r + /T_ = 75 (3^7) as a result of their 



Monte Carlo analysis (this is the most recent result known to us). 

Since we expect the accuracy of our results at q = 3,4 to be comparable with that 
found at q = 2, the main source of uncertainty for the extrapolated results is in the 
extrapolation itself. Comparison between the extrapolated and the Monte Carlo values for 
£oV£o" shows that our error does not exceed 10% in this case. We find quite reasonable to 
assume the same level of accuracy for our prediction on r + /T_. We summarise in Table 5 
the situation about the amplitude ratios in two-dimensional percolation we considered in 
this paper. 



12 The existing e-expansion results are unreliable in two-dimensions since the upper critical dimension 
in the problem is d c = 6. 

13 The value A + /A- = —1 is quoted in Refs. [B8| and 0. We do not understand the origin of this 
discrepancy. 

14 It is immediately checked that extrapolating in q gives essentially the same results for £q/£q and 
r + /r_ while the value of Rc gets modified by less than 10%. 

15 Starting apparently from the same lattice amplitudes, Aharony |38| finds Rc ~ 4.1 — 4.2 and this 
result is quoted also in |7J]. 
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6 Conclusion 

In this paper we used the form factor bootstrap approach to compute several universal 
quantities for the g-state Potts model and isotropic percolation in two dimensions. The 
results have been obtained truncating the spectral series for the two-point correlators 
at the level of the two-kink contribution. We showed by comparison with exact results 
that this approximation is sufficient to provide remarkably accurate values for physical 
quantities like central charge, scaling dimensions and critical amplitudes. In particular, 
we gave the first theoretical prediction of the universal amplitude ratios in the g-state 
Potts model for q = 3, 4. 

Our ability to provide precise predictions for percolation is limited by the difficulty to 
solve Eq. ( ^.42[ ) around q — 1. We gave accurate results for some of the amplitude ratios 
and relied for the others on an extrapolation based on the values obtained at q — 2, 3, 4. 
Our predictions (including those from extrapolation) are found to be in good agreement 
with the existing lattice estimates when the latter are able to provide reasonably accurate 
results. This is not the case for the ratio of the mean cluster size amplitudes above and 
below the percolation threshold, for which the different lattice determinations span more 
than one order of magnitude. In light of this, our extrapolated value for this quantity 
represents substantial progress. 

All the results of this paper have been obtained within the form factor framework, 
using the S'-matrix as the only input and without relying on data coming from other 
approaches (e.g. TBA). This was possible also because we computed universal combi- 
nations of amplitudes concerning two renormalisation group trajectories related by dual 
symmetry, and then describable through the same scattering theory. Other universal 
ratios can be defined which involve also amplitudes computed at the critical tempera- 
ture but in presence of an external magnetic field. When trying to determine them in 
integrable field theory, one faces the problem of fixing the relative normalisations of am- 
plitudes computed through different integrable scattering theories. It is remarkable that 
also this problem can be solved using the results of the TBA |42|, ^] and some more 



recent developments P5 , |4"4| ] . It was shown in Ref . [ S3 1 through the basic example of the 
Ising model how the universal ratios at nonzero magnetic field can be computed along 
these lines. 
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Appendix A 

Consider the matrix element 

(K aP (e')\m\Kf3a(0)) • (A.l) 

It can be related by analytic continuation to a matrix element between the vacuum 
and a two-kink state. There are, however, two different ways to perform the analytic 
continuation, corresponding to the fact that the final two-kink state can be an "In" or 
"Out" asymptotic state. This leads to the two equations 

(K aP (e')mo)\K 0a (e)) = (O a mO)\K a0 (8' + i7r)K Pa (e)) + 2n5(6 - 9')(0 a \®\0 a ) , (A.2) 

(K a0 (e')\$(O)\Kp a (O)) = (O mO)\K 0a (9)K a0 (9' - in)) + 2n5(9 - 9')(Op\®\O ) , (A.3) 

where the delta function terms are disconnected parts which take into account kink- 
antikink annihilation. While these equations are kinematically identical to those one 
obtains when dealing with ordinary particles, the important difference to be noticed is 
the appearance of two different vacuum states. Subtracting the second equation from the 
first, one gets 

(O a mo)\K a p(9' + i7r)K Pa (9)) = (0^(0) \K 0a (9)K a0 (9' - m)) 

+ 2n6(9 - 0')«O^|$|O^) - (0 Q |$|0 Q » . (A.4) 



As long as 9 ^ 9', this is exactly the crossing relation ( |3.3| ). For 9 = 9', the above 
equation implies the presence of a simple pole in the two-kink form factor with residue 

Res e=e/ (0 a m0)\K al3 (9' + iix)K Pa {9)) = -K^eMM^)\Kpa{9)K aP {9' - in)) 
= i((0 a \$\0 a ) - (0f,\$\0p)) . (A.5) 

The presence of a pole in the two-particle form factors when the rapidity difference equals 
m is known to reflect the non-locality of the operator with respect to the fields which 
interpolate the asymptotic particles (see e.g. |T3j]). We see that in the low temperature 



formalism this amounts to the fact that the operator has different expectation values 
on different vacua. Among the operators considered in this paper, this is the case of 
the magnetisation cr a (x) which is non-local with respect to the disorder operators which 
create the kinks. 
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Table Caption 



Table 1 . Central charge and scaling dimension of the magnetisation operator in the 
g-state Potts model. The results obtained through the form factor approach in the 
two-kink approximation are shown below the exact values. 

Table 2 . Values of some amplitude combinations in the g-state Potts model as obtained 
through the form factor approach in the two-kink approximation. 

Table 3 . Universal amplitude ratios in the g-state Potts model. The exact result for 
A + /A- follows from duality; the other values are computed through the form factor 
approach in the two-kink approximation. 

Table 4 . Series and Monte Carlo estimates of amplitudes for bond percolation on the 
square lattice and site percolation on the triangular lattice. The results marked by 
the superscript a, b and c are taken from Refs. ||37|| , [f40l and pOU , respectively. 



Table 5 . Universal amplitude ratios in two-dimensional percolation. The dagger signals 
extrapolated values. The results marked by the superscript a, d and c are taken 
from Refs. |37j , |7j and [^] , respectively. The series/MC estimates for R c and Rt 
are obtained using the results of Table 4. 
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Figure Caption 



Figure 1 . Schematic representation of the two-body scattering amplitudes in the Potts 
model scattering theory. The continuous lines represent a kink, the dotted lines the 
bound state B. 

Figure 2 . The three-particle vertices in the Potts model scattering theory. 

Figure 3 . Schematic representation of a two-kink form factor. 

Figure 4 . (a) Two-kink form factor of the energy operator at the resonant rapidity 
difference 9 = 2m; (b) two-kink form factor of the disorder operator at the resonant 
rapidity difference 9 = 2in/3. 

Figure 5 . Central charge in the g-state Potts model. The continuous line is the exact 
formula (|2.6| ) . For q < 3 the dotted line gives the result of the c-theorem sum rule 
computed through the form factor approach in the two-kink approximation. For 
q > 3 the bound state B enters the physical spectrum. The upper dotted branch 
represents the sum of the kink-kink contribution (lower dotted branch) and the 
single particle bound state contribution. 
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4 


c 


1/2 

1/2 


4/5 
0.792 


1 
0.985 


X(j 


1/8 

1/8 


2/15 
0.128 


1/8 
0.117 



Table 1 



Q 


1 


2 


3 


4 


T + ml/B 2 
£o~™o 


1.001 


2 
1 


0.973 
0.998 


0.476 
0.992 
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Q 


2 


3 


4 


A + /A_ 


1 


1 


1 


r+/r_ 


37.699 


13.848 


4.013 


Re 


0.3183 


0.1041 


0.0204 


so /so 


3.162 


2.657 


1.935 


Rt 


0.3989 


0.3262 


0.2052 
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SQ bond 


TR site 


K 


-4.24° 


-4.37° 


B 


0.77 6 


0.78 6 


f+ 


0.134 6 


0.128 6 


u 


- 


0.313 c 
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This work 


Series/MC 


A + /A_ 


1 


l a 


r+/r_ 


74.2+ 


14 - 220 d 


SO /SO 


3.76+ 


4.0±0.5 C 


Re 


3.0+ 


2.7-2.8 


Rt 


0.926 


1.1 



Table 5 
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